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ABSTRACT 

We present spectra of 1142 colour-selected stars in the direction of the Sagittarius 
Dwarf Spheroidal (Sgr dSph) galaxy, of which 1058 were taken with VLT/FLAMES 
multi-object spectrograph and 84 were taken with the SAAO Radcliffe 1.9-m telescope 
grating spectrograph. Spectroscopic membership is confirmed (at >99% confidence) 
for 592 stars on the basis of their radial velocity, and spectral types are given. Very 
slow rotation is marginally detected around the galaxy's major axis. We identify five S 
stars and 23 carbon stars, of which all but four carbon stars are newly-determined and 
all but one (PQ Sgr) are likely Sgr dSph members. We examine the onset of carbon- 
richness in this metal-poor galaxy in the context of stellar models. We compare the 
stellar death rate (one star per 1000-1700 years) to known planetary nebula dynamical 
ages and find that the bulk population produce the observed (carbon-rich) planetary 
nebulae. We compute average lifetimes of S and carbon stars as 60-250 and 130-500 
kyr, compared to a total thermal-pulsing asymptotic giant branch lifetime of 530-1330 
kyr. We conclude by discussing the return of carbon-rich material to the ISM. 

Key words: stars: AGB and post-AGB — stars: carbon — stars: mass-loss - 
circumstellar matter — stars: abundances — galaxies: individual: Sgr dSph 



1 INTRODUCTION 

The mass return from stars to the interstellar medium (ISM) 
is a powerful driver of chemical change in a galaxy. The 
abundance and amount of material returned to the ISM de- 
termines the chemical make-up of successive generations of 
stars and planets. A large fraction of the metals returned 
to the ISM are ejected by highly-evolved asymptotic gi- 
ant b ranch (AGB) stars (e.g. iGehr j[l9 89: Zhukovsk a et al.l 
2008) which lose mass via pulsation-enhanced, dust-driven 
winds (e.g. IWoitkd l2006aft . Dusty material returned by 
oxygen-rich stars include silicates, primarily pyroxenes and 
olivines (Mgi-^Fe^SiOs, Mg2-2a;Fe2iSi04), oxides such 
as AI2O3, and p o tentially metals and metal alloys (e.g. 
Sloan et all 120031; Ivan Loon et all 120061: IMcDonald et all 



20091: iTamanai et al.1 120091 : IMcDonald et all l20ld . l2011bl : 



Norris et all l2012t ). Carbon stars return primarily amor 



phous carbon, silicon carbide, and other car b on-rich species 
|Treffers fc Cohenl 1 1974 iMillar et al.1 l2000l : lLagadec et ail 
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12007ft . Carbon-rich dust species tend to have signifi- 
cantly higher opacities than oxygen-rich species, mean- 
ing that stellar outflows can be driven more effectively 
l|Woitkel [2006b). It has been proposed that the 'super- 
wind' (the final, high-mass-loss-rate phase of the AGB 
wind that ejects the remainder of the star's atmo s phere 
and ends its nuclear-burning life; van Loon et al.1 1 19991 : 
IWillsonl l2000l : iFerrarotti fe Gail 12006ft may sometimes be 
triggered by the increase in C/O ratio (|Winters et al.ll2000l : 
lLagadec fc ZiUstrdl200Sft . 

Giant stars undergo dredge up of fusion-processed ma- 
terial from their cores at key stages in their post-main- 
sequence evolution (e.g. ISmith fc Lambert! Il985l : iMowlavil 
1999; iHerwid | 20041 ') . Carbon is an important element in 
the dredge-up process, particularly during third dredge- 
up, which occurs on the thermally-pulsing AGB (TP- AGB) 
|Drav et alj|2003ft . It is over-abundant in the dredged mate- 
rial, which can result in an increase in the carbon-to-oxygen 
(C/O) ratio from a main-sequence value of C/O « 0.4 to 
C/O > 1. Oxygen and carbon are among the most abun- 
dant metals, and CO is one of the first molecules to form. 
This binds away free carbon and oxygen, leaving the more 
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abundant element to dominate the stellar chemistry. When 
C/O exceeds unity (by number), the star transforms from 
an oxygen-rich star to a carbon-rich star. The length of this 
carbon-rich period (which lasts for the remainder of the 
star's nuclear-burning life) and the mass loss that occurs 
during it, directly affects how much carbon is returned to 
the interstellar medium. 

Of particular importance for the historical enrichment 
of the ISM, in our own and other galaxies, is the varia- 
tion of the onset and length of this carbon-rich phase with 
metallicity and initial stellar mass. Only some stars become 
carbon-rich: low-mass stars do not dredge up sufficient car- 
bon to become carbon-rich, while high-mass stars undergo 
hot-bottom burning via the CNO cycle, which reduces the 
carbon abundance in the dredged-up layers , limiting surface 
carbon enrichment (lLattanzio et al.lll996l ). Exactly where 
these limits lie depends on t he complex parameters of stel- 
lar mass loss and dredge-up (| Ventura fc Marigoll2010l ). For 
stars of (initially) solar abundance, the bounds for becom- 
ing c arbon rich lie at initial m asses of approximately 2 and 4 
Mq l|Marigo fc Girardill2007l ). Lower-metallicity stars, which 
have lower initial oxygen abundances, require less carbon 
to be dredged up to become carbon-rich. The mass range 
of low-metallicit y carbon stars is therefore expected to be 
much larger (e.g. iMatsuura et a.1 .1120051 : IZiilstra et al.| [2006b ; 
iGirardi fc Marigoll2007l )~ 

The onset of carbon-richness at different metallicities 
and masses can be probed by examining stellar populations 
with different metallicities and ages. Among the most in- 
teresting is the Sagittarius Dwar f Spheroidal (Sgr dSph) 
or E lliptical Galaxy (Sag DEG; llbata et all 1 1994 1 19951 . 
1997): this metal- poor galaxy l i es at a relatively-clo s e dis- 
tance of 25 kpc jMateo et al.l 1 19951 : iMonaco et"afl |2004 
iKunder fc Chaboverll2009T)" and is being tidally disrupted by 
the Milky Way jTotten fc Irwirj|l998l ; iMaiewski et al.ll20031 ; 
lLaw et al]|2005h . The galaxy contains several identified pop- 
ulations. While the delineation between these populations is 
indistinct, at least at present, they have been broadly sum- 
marised bv lSiegel et al.l (|2007T l into five separate groups: 



(i) a very metal-poor population with [Fe/H] = -1.7 dex, 
[a/Fe] w +0.2 dex, and an age of t « 13 Gyr, which derives 
from (and is mostly contained within) the globular cluster 
M54; 

(ii) a sizable, metal-poor population with [Fe/H] as -1.2 
dex, t = 10-12 Gyr, most p revalent in the outskirts (see also 
lLavden fc Saraiedinill2000l ); 

(iii) a bulk population with a spread of metallicity be- 
tween [Fe/H] = -0.7 to -0.4 dex, [a/Fe] « -0.2 dex, and 
ages between t = 4 - 8 Gyr, which d ominates in the galaxy's 
core (see also iBellazzini et al.l 120061 ); 

(iv) a small, higher-metallicity population ([Fe/H] « 



-0.4 to -0.1, 



t/F 



+0.2 dex, t 



Gyr; Saraiedini fc Lavdenlll995l; l Lavden fc Saraiedinill2000l ; 



IZiilstra et alj|2006al : ISiegel et al.ll2007f ); and 

(v) a probable, very small, metal-rich population ([Fe/H] 
~ +0.5, [a/Fel ~ dex, t < 1 Gyr, newly discovered by 
ISiegel et alj|2007h . 



Several globular clusters may also be historically associated 
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Figure 1. 2MASS colour-magnitude diagram of the Sgr dSph 
core, showing the colour cuts used to isolate the Sgr dSph RGB 
(region (a)) and AGB (region (b)). Region (c) is inhabited by 
dusty AGB stars. Region (d) contains foreground dwarf stars 
and background galaxies. Region (e) contains a population of un- 
known providence, which is discussed in Section 33] The dashed 
blue lines show the regions (f,g) used to calculate RGB and AGB 
star counts. 



with the galaxy and its tidal strearrQ; these have metallic- 
ities ibet^enJFe/H] = -2.16 (Terzan 8) and -0.32 (Terzan 
7) (|Harrij[i99^ ). The best-studied of these is M54, which 
lies at the projected centre of the galaxy. 

The Sgr dSph's most-evolved AGB stars should be close 
to the lower mass limit for carbon-richness. Previous works 
have noted several carbon stars i n the galaxy, and linked 
them to the metal-ri ch population (La gadec fc Ziilstrall2008l ; 
lLagadec et~aill2009h . In this work, we examine spectra of a 
significant fraction of the galaxy's giant stars, in order to 
confirm their membership of the galaxy and determine the 
point at which stars become carbon rich. 



2 OBSERVATIONS 

The observed stars were sampled from the 2MASS catalogue 
(|Cutri et al.l 120031 ). The 2MASS (J - K s ) vs. K 3 colour- 
magnitude diagram (CMD) for the 3° around M54 (which 
we take to be the centre of the Sgr dSph) shows a dis- 
tinct feature associab le with the Sgr dSph giant branch 
ijMaiewski et al ] |2003l : see also Figure [1} . Stars represented 
by this feature can be separated from the foreground Galac- 
tic population using a colour cut of K a > 20.35 — 9( J — K B ) 
mag (left-most diagonal line, Figure [T]), which works toler- 
ably well for stars of K s < 12 mag. Priority was given to 
stars above the RGB tip (K s ~ 10.7 mag). 

The limits of allowed fibre positions and field of view of 
the FLAMES spectrograph meant that these observations 
contained some stars above the colour cut. These are likely 



1 Though authors differ on which clusters are associated, the 
following are generally accepted: M54, Terzan 7, Terzan 8, 
Arp 2, Palomar 2 (still uncertain), Palomar 12, W hiting 1 
ijda Costa fc Armandrofflll995l : I Law fc Maiewskll2010l and refer- 
ences therein). 
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Table 1. Details of observations taken. 



Telescope 



Field Centre (J2000.0) 
RA Dec 



Observation Instrument Coverage Resolution Number of 
date (nm) (A/AA) objects 



VLT 



S A AO- 1.9m 



18 56 24 
18 56 27 
18 59 39 
18 54 00 
18 53 24 
18 58 12 
18 55 12 
18 48 29 
18 52 12 



-30 31 48 
-31 22 47 
-31 01 10 
-30 44 59 
-30 20 59 
-30 42 00 
-30 12 00 
-30 12 02 
-29 59 58 



2010- 
2010- 
2010- 
2010- 
2010- 
2010- 
2010- 
2010- 
2010- 



06-24 
09-20 
09-20 
09-22 
09-23 
09-24 
09-25 
09-26 
09-27 



2009-07-04 to 

2009- 07-08 

2010- 07-28 to 
2010-08-04 



GIRAFFE 



UVES 



SpCCD 
SpCCD 



574.1-652.4 
643.8-718.4 
638.3-662.6 
693.7-725.0 
476.0-684.0 



400-730 



7400 
8 600 
28 800 
23 900 
47 000 



1 750-3 170 



642-722 15 000-16 800 



1015 



43 



28 
262 
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Figure 2. Spatial distribution of the observed VLT/FLAMES 
fields, against stars from the 2MASS catalogue with (J — K B ) 
colours matching the Sgr dSph giant branch. The open circles 
show the SAAO targets. The large magenta cross marks the glob- 
ular cluster M54, the galaxy's nominal centre. Red points denote 
radial velocity members, blue points denote radial velocity non- 
members, green points denote stars where a radial velocity was 
not accurately determined. 

Galactic stars, which are useful for determining the Galac- 
tic velocity spread in the region, but these stars may also 
include some unusual objects from the Sgr dSph. We show 
the spatial location of our observations within the Sgr dSph 
in Figure [2] noting that they are very much concentrated in 
the galaxy's core (Figures 

The bulk of our observed spectra were obtained us- 
ing the ESO VLT/FLAMES+MEDUSA multi-fibre spectro- 
graph in GIRAFFE+UVES mode. These cover nine fields 
of the Sgr dSph core, each taken in four settings (L714.2, 
L682.2, H651.5B and H71O.50- Simultaneous spectra were 
taken using the UVES instrument using the red arm, cross- 
disperser CD3 and a central wavelength of 580 nm. The 
higher resolution of UVES allowed us to focus on objects 
near the AGB tip. We note, however, that the echelle for- 



2 Full data on each setting is available in the FLAMES manual: 
http: / / www. cso. org/sci / facilities / paranal / instruments / flames / 
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Figure 3. Location of colour-selected RGB stars (grey dots, re- 
gion (a) in Figure [TJ. Overlaid as small blue points are (near-) 
naked AGB stars, selected from region (b) of Figure [T] Overlaid 
as larger red points are dusty AGB stars, selected from regions 
(c) and (e) of Figure [T] The blue, contaminating population to 
the top right is the Galactic Disc, the group at RA = 285° to 
290°, Dec = -37° is the Galactic star- forming region associated 
with NGC 6727. The black lines mark the regions we use to define 
the galaxy's "core" and "main body" within the text. 



mat of UVES and the faint nature of many of our targets 
means that the spectra suffer from substantial distortion 
where spectral orders meet, limiting our ability to measure 
the shape of broad molecular features. 

The large size of the Sgr dSph compared to the 
FLAMES field of view (Figures \3§ prompted observations 
of additional stars using the grating spectrograph on the 
SAAO Radcliffe 1.9-m telescope over two runs. Observations 
were taken during variable conditions, however, leading to a 
variable signal-to-noise in the final data. TiO or CN bands 
could be discerned in 69 objects. 

These observations are summarised in Table [1] 
The VLT data were extracted from the standard 
pipeline reduction provided by ESO. As the GIRAFFE data 
consist of single exposures in each setup, it was not pos- 
sible to mitigate against cosmic rays. The data were sky 
subtracted, but not corrected for telluric absorption as this 
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is not required for our current analysis. The simultaneous 
UVES data consist of four exposures per setup. The four 
resulting spectra were median combined to remove cosmic 
rays. 

The SAAO data were reduced following the normal pat- 
tern of bias and dark subtraction. Flatfielding was not per- 
formed due to the poor quality of fiatfields obtained during 
the observing run. Differences in the responses of adjacent 
pixels in the detector are small compared to the noise in 
the final spectra; while lack of flatfielding means we cannot 
determine the overall shape of the spectra, it should not sig- 
nificantly affect the accuracy of continuum-divided spectra. 
Literature data for 324 objects were als o sourced 



from Ibata et alj 



Bonifacio et alj 



19951); IWhitelock et all (1 199311 
( 2001 12001); iLanfranchi et all d2Q06h 
McDonald fc van Loonl (f2007h: lLagadec et all (|2009T ) 
Giuffrida et all d,201Ch : ICarretta et al.1 ll2010l)~lhe Genera l 
Catalogue of Galactic Carbon Stars (|Alksnis et al.l 200ll); 
the General Catalogue of Variable Stars llSamus et al. | 
2006); and the Radial Velocity Experiment (|Siebert et al.l 
20111 ). These sources contain a variety of radial velocity and 
spectral typing data which we use in later figures. 



3 ANALYSIS 

3.1 Spectral typing 

An initial visual inspection of each spectrum was performed 
to separate carbon- and oxygen-rich stars for the purpose 
of radial velocity determination. This was done on the basis 
of C2/CN and TiO/VO bands, respectively (see Figure |4}. 
The final spectral types are listed in Tables [5] and [3] They 
are displayed along with literature results in Figure [5] 

VLT spectra without discernable molecular bands were 
assumed to be of spectral type K. SAAO spectra with no dis- 
cernable bands were labelled as unknown spectral type (U). 
The difference in approach arises from the different data 
qualities. Due to poor observing conditions during the 2010 
observing run, some SAAO data are of such low signal-to- 
noise in some cases that even significant molecular bands 
would not be visible. The VLT spectra are of sufficient 
signal-to-noise that TiO bands would always be visible if 
they are present. The C2 and CN features of carbon stars 
are not as clear as TiO features in the VLT data, meaning 
some carbon stars could theoretically be misclassified as K 
stars. However, the stars we classify as K stars are all at 
low luminosities where only a few extrinsic carbon stars are 
expected, so we believe we have listed all observed intrinsic 
carbon stars in our sample. 

One galaxy (#29) is observed in the sample. The spec- 
trum appears featureless, but optical and 2MASS images of 
the source show that it is elongated. We therefore discount 
it from further analysis. This leaves us with 718 K-type, 
290 M-type and 7 C-type stars observed with GIRAFFE, 
respectively; 1, 41 and 1 with UVES; and <15, 54 and 15 
with the SAAO 1.9-m. The higher fraction of carbon stars 



3 We do not include Whitelock et al.'s WMIF 9 and 16-27 as, 
although they are radial velocity members of the Sgr dSph, they 
are out of our spatial region of analysis. Additionally, WMIF 13 
cannot be identified in the 2MASS catalogue. 



in the SAAO data is an intentional bias, introduced by pref- 
erentially observing the brightest AGB stars. 



3.2 Radial velocity determination 

Radial velocities for the VLT/GIRAFFE targets were de- 
termined using the Starlink routine HCROSS and corrected 
from topocentric to barycentric co-ordinates using the task 
RV. We used the stars 1013, 1014 and 1015 were used as tem- 
plates for the K, M and C stars, respectively (Figure[4]|. The 
three templates were visually cross-correlated with a 4000 K 
MARCS model with log q = 1.5 dex, [Z/H] = -0.75 and [a/FeJ 
= +0.3 (|Gustafsson et alJl97rJ . l200alMcDonald et alj|2009h 
to provide absolute radial velocities with an estimated error 
of 2 km s" 1 . 

Each star was cross-correlated against its relevant tem- 
plate over 30 spectral regions across the four spectral set- 
tings. These spectral regions were chosen to bracket par- 
ticular features whilst avoiding areas of telluric absorption. 
Values lying more than 2.17cr (97% confidence interval) from 
the median were iteratively clipped, and a single radial ve- 
locity for each star was determined using the median of these 
remainder. The resulting random error in radial velocity av- 
erages 1.8 km s" 1 , with 87% of sources having errors <10 
km s _1 . The remaining 13% of sources are mostly of too low 
signal-to-noise to accurately determine a radial velocity. 

Lacking an absolute radial velocity standard, we assume 
that the average Sgr dSph star (those with 100 < v T < 200 
kms -1 and ov < 10 km s" 1 ) has a radial velocity of v r = 141 
km s" 1 (cf.llbata et al.lll997l; Izilistra fc Walshll 19961 ; iHarris] 
ll996l ; lGiuffrida et alj|20ld : iPenarrubia et al.ll201lD . and off- 
set our velocities accordingly. Only one of our stars has a 
published rad ial velocity (2MASS 18554672-3035248, our 
source #346; ICarretta et~al"1 12010| . their source 23001277, 
v t = 147.6 km s" 1 ). Our velocity (148.22 ± 1.05 km s" 1 ) is 
in agreement within the errors, showing that our absolute 
velocity calibration is accurate. 

The radial velocities of the 43 observed UVES targets 
were determined using a similar method. Here, however, ra- 
dial velocities could not be calculated using the Starlink 
HCROSS package. This appears to be due to problems caused 
to the Fourier transformed data by poor order matching, in- 
creased noise in the spectrum, and the generally richer and 
more- varied spectra of these typically more-evolved stars. To 
combat this, we removed small-scale structure in the spec- 
tra by smoothing them by a Gaussian of a = 10 pixels and 
removed large-scale structure by dividing them by the same 
spectrum smoothed by a Gaussian of a = 100 pixels). Each 
spectrum was sliced into seven segments, each of 3000 pixels, 
and each segment was cross-correlated with the correspond- 
ing segment from every other star. As above, an iteratively- 
clipped (2a) median was used to calculate a unique radial 
velocity and the standard deviation of the measurements 
used to create a measurement of the velocity error. 

Of the resulting 43 radial velocities, 42 are very tightly 
grouped and obviously belong to the Sgr dSph. The 43rd 
star, #1033, is Galactic. We adjust the radial velocity zero 
point so that the average of the 42 Sgr dSph stars lies at 141 
km s _1 to place them on a common scale with the GIRAFFE 
spectra. 

Accurate radial velocities could not be determined for 
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Figure 4. Examples of typical spectra, selected from stars just below the RGB tip (top: M-type #105, middle: C-type #1015, bottom: 
K-type #949), in the four observed filter settings. Notable atomic and molecular lines are shown. 

Table 2. Spectral types and radial velocities of targets observed with the VLT. Stars with IDs ^1016 were observed 
with UVES, the remainder with GIRAFFE. The right-hand five columns denote the median heliocentric radial velocity 
determined from the spectra, the standard deviation of radial velocities determined for that object (which we take as the 
error), the probability of membership of the Sgr dSph, and the TiO and ZrO band strengths as described in the text. A full 
table is available in the online version. 



ID 


RA 


Dec 


Spectral 


2MASS 




"helio 




P 


TiO 


ZrO 




(J2000.0) 


(J2000.0) 


Type 


K a {J - 


K s ) 


(km s" 1 ) 


(km s _1 ) 


% 


index 


index 


0001 


18 55 12.44 


-30 23 39.0 


K 


9.92 


0.35 


29.53 


12.19 


0.00 


1.04 


0.91 


0002 


18 55 18.97 


-30 21 01.7 


K 


9.26 


0.75 


-14.93 


1.24 


0.00 


1.03 


0.95 


0003 


18 54 12.34 


-30 24 10.0 


K 


12.28 


0.78 


-8.96 


1.68 


0.00 


1.01 


0.96 



the SAAO data due to the low spectral resolution and poor 
data quality. 



3.3 Membership probability 

To define a membership probability of the Sgr dSph, we ap- 
proximate the 556 stars with Sgr-dSph-like velocities (as de- 
fined above) and 312 Galactic foreground stars (those with 
— 150 < v T < 100 km s _1 and a v < 10 km s" 1 ) as two 
Gaussians, with < vs g r >= 141.0 and o"s gr = 12.7 km s _1 , 



and <vcsi >= —10.3 and o Q a i = 45.3 km s _1 , respectively. 
Membership probability of the Sgr dSph (P) was then de- 
fined as: 



P = 



556exp(- (, '- < g > >)2 ) 
v 312exp^- ( "-<^' >)2 ^ 



■ exp 



{v 



Av 



Av 



,(1) 



© 9999 RAS, MNRAS 000.HTfT4l 



6 I. McDonald et al. 







1 2 3 

(J-K s ) 


4 




1 


I . . I I I I 

X 


1 


X" 












■ 




X 


1 


I I I I I 


1 





9.5 
10 

10.5 
11 

11.5 



0.85 0.9 0.95 



1 1.05 1.1 1.15 1.2 1.25 
ZrO 



»G 
G 



• -6 



______ *AP>_ _ . 



0.8 0.9 



1.3 1.4 




1.5 1.6 



Figure 5. Top panels: Colour-magnitude diagrams showing the distribution of spectral types in the observed sample (circles) and 
literature (plus signs). Colours represent: (red) carbon stars, (green) M giants, (blue) earlier spectral types, (yellow) undetermined or 
unlisted spectral type. Hollow symbols show SAAO targets, smaller filled symbols show FLAMES targets, larger symbols show UVES 
targets. The observed background galaxy and Galactic carbon star PQ Sgr are labelled. The right-hand panel shows a zoomed-in region 
of the left-hand panel near the RGB tip, which is marked by the horizontal line. 

Bottom panels: Molecular band strengths of the GIRAFFE sample, using indices from Otto ct al. (2011J). Carbon stars are shown by 
circles, the S stars by crosses. Bottom-left panel: The colour scale shows (J — K B ) colour, with darker stars being redder. Bottom-right 
panel: Colour-magnitude diagram, showing the colour cuts applied from Figure [T] The colour scale shows the strength of the ZrO band 
index. Carbon and S-type stars in the UVES and SAAO samples are shown as lighter (red) points. 



Table 3. Spectral types of targets observed with the SAAO Rad- 
cliffe 1.9-m. A full table is available in the online version. 



ID 


RA 


Dec 


Spectral 




2MASS 








Type 


K s 


(J - K s ) 


2001 


18 56 11.24 


-29 03 54.9 


C 


06.31 


1.92 


2002 


18 53 59.53 


-29 58 57.4 


C 


09.26 


1.52 


2003 


18 54 09.02 


-30 23 20.2 


M 


09.28 


1.29 



where v is a test radial velocity, v T is the calculated (clipped 
median) radial velocity and a v is its error. 

Figure [6] shows a histogram of these velocities. A dis- 
cernable break is apparent between Galactic stars and Sgr 
dSph members at Whoiio ~ +100 km s _1 , beyond which most 
stars are Sgr dSph members. The velocity selected (non- 
)members and stars whose membership could not be reliably 
determined are shown in Figure [7] which shows that the Sgr 
dSph members clearly lie almost exclusively along the red- 
der giant branch. Based on our probability analysis, we find 



593 stars are members of the Sgr dSph at high (>95%) con- 
fidence, 338 stars are members of our Galaxy, and 127 stars 
could not reliably be placed in either category. 

3.4 Band indices and S stars 

Having a spectrum in the stellar rest frame, we can now per- 
form more accurate spectral a nalysis. To do this , we adopt 
the TiO and ZrO indices from lotto et al.l (|201ll ). inverting 
the TiO index as per their figure 2. Namely, these are the ra- 
tio of the average flux at 640-646 nm to that at 647.5-653.5 
nm for the ZrO band, and 696.5-702.8 nm to 706.5-717.5 
nm for the TiO band. We do not cover the spectra l regio n 
required to calculate the CT51 index of lotto et alj (|201ll ). 
so we cannot repeat their analysis exactly. 

Figure[5]shows the band indices we derive. Objects with 
IztO > 1-05 were visually searched for ZrO bands. We iden- 
tified S stars by com parison to carbon-enriched spectra from 
Ivan Eck et alj (|201ll ) . As these spectra are not yet fully pub- 
lished, a visual comparison to their Figure 2 was necessary. 
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Figure 6. Histogram of derived radial velocities of FLAMES tar- 
gets. The dashed curves are Gaussians fit to the Sgr dSph and 
Galactic populations, as described in i|3.3l The Galactic popula- 
tion is relatively well separated, lying almost exclusively below 
100 km s —1 (vertical line). 
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Figure 7. Colour-magnitude diagram showing the positions of 
Sgr dSph members (red), Galactic foreground stars (blue) and 
stars of indeterminate location (green) observed with the VLT, 
overlaid on a density plot of 2MASS sources on a logarith- 
mic greyscale. Circled objects highlight stars with velocities co- 
incident with the Sgr dSph, but colours that are not. These may 
be metal-poor or post-AGB member stars. 

Four M stars were reclassified as S stars identified on this 
basis: #361, 424, 909 and 942. A fifth star from the UVES 
spectra (#1054) was also reclassified as an S star. We return 
to these stars in Section \4. 51 

The remainder of the stars cluster close to unity, with 
the cooler, M-type, stars displaying increasingly large values 
of the TiO index. The carbon stars, although equally cool, 
do not share the TiO or ZrO features, thus are the only very 
cool stars to appear near unity in both indices. 



4 DISCUSSION 

4.1 Velocity gradients 

To investigate radial velocity gradients in the Sgr dSph, we 
have taken the K star population in each observed field and 



computed an average velocity for those K stars with veloc- 
ity uncertainties under 10 km s _1 and corrected these back 
to the Galactic standard of rest using the NASA/IPAC Ex- 
tragalactic Database (NED) velocity converteiQ We show 
these in Figure [8] A possible north-east-south- west gradient 
of ± ~2 km s _1 exists, which could conceivably be aligned 
with the galaxy's minor axis and/or parallel to the Galactic 
plane, however it is at the limit of detection and therefore 
cannot be conclusively shown. Neither can we determine on 
the basis of these data alone whether such a gradient is in- 
dicative of a tidal distortion or a true rotation. 

4.2 Investigating the metal-poor population 

Figure [7] shows nine objects with a high membership prob- 
ability (>99%), but with (J — K a ) colours much bluer than 
the main Sgr dSph population, namely #10, #66, #105, 
#155, #170, #190, #232, #243 and #260. Given the slight 
overlap between the Galactic and Sgr dSph velocity distri- 
butions (Figure [6}, we cannot be certain that these are Sgr 
dSph stars based on their velocities alone. These stars are 
important, as they may be members of the metal-poor pop- 
ulation, or post-AGB stars. 

Following a thorough visual investigation of their spec- 
tra and associated literatue data, we are unable to conclu- 
sively determine (non-)memberships of these stars. However, 
we suggest that #155 and #260 are probably Sgr dSph 
stars, based on their proximity to the Sgr dSph branch on 
the (J-K a ) colour-magnitude diagram, and on the carbon- 
richness of #260. We further suggest that #66, #105 and 
probably #243 are Galactic, based on their proper motions 
IjRoeser et al.ll2010l ). The other stars remain at an undeter- 
mined distance, as we have insufficient data to determine 
their location without a full spectral synthesis to determine 
their chemical abundances. Removing #66, #105 and #243 
from the Sgr dSph sample has little effect on our overall re- 
sults (entirely within the measurement errors). The most no- 

4 http:/ /ned. ipac.caltech.edu/help/velc_help. html 
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table change is that removing #243 from the most-northerly 
field decreases its average heliocentric radial velocity by 0.9 
km s _1 and its standard deviation by 0.4 km s" 1 (cf. Figure 



4.3 The stellar evolution rate 

4-3.1 Determining evolutionary rates 

One can, to first order, estimate the rate at which stars 
evolve off the AGB by comparing the observed colour 
magnitude diagram to stellar evolutionary tracks. RGB evo- 
lutionary models tend to be more easily-available and much 
better calibrated than those of the AGB. Furthermore, many 
evolutionary models (including the models used here) do 
not include mass loss, or assume a prescription of mass loss 
which may not be appropriate. For stars above ~1 Mq, RGB 
mass loss is relatively small and can usually be ignored, and 
(as we will show) RGB evolutionary rates have little depen- 
dence on stellar mass in most normal regimes. It is therefore 
easier to use RGB models to determine the RGB evolution 
rate, then determine the AGB evolution rate by using the 
ratio of star counts on the RGB and AGB. This should be 
identical provided every RGB star eventually ends its AGB 
evolution above the RGB tip. 

Throughout this process, we will refer to the regions la- 
belled (a) through (g) in Figure [T] Star counts for these re- 
gions are displayed in Table [4] In this Figure, region (a) cor- 
responds to the unreddened RGB, region (b) to the unred- 
dened AGB and region (c) to the reddened AGB. Region 
(d) is largely devoid of interesting stars, containing mostly 
background galaxies. Region (e) may correspond to a red- 
dened RGB, though the true identity of these stars (be they 
RGB or AGB, carbon- or oxygen-rich) is unknown. Regions 
(f) and (g) show the RGB and AGB populations we use to 
determine evolutionary rates. Note that these regions, each 
covering AK S = 0.4 mag, avoid the 0.1 A^-band magnitudes 
either side of the RGB tip (at A" s = 10.68 mag) to account 
for uncertainties in the 2MASS photometry and uncertainty 
in the position of the RGB tip itself. 

Table [4] also contains four columns which extrapolate 
the number of carbon stars and Sgr dSph members to the 
entire 3° x 3° core region (as defined in Figure [31). In each 
case, we show two numbers. The first is the 68% confidence 
interval for the number in that region, based on multiply- 
ing the relevant observed fraction (carbon star fraction or 
membership fraction) by the raw counts. The second shows 
the same confidence interval taking into account the Pois- 
sonian error in the number of stars. This latter correction 
is necessary when exploring evolutionary rates, as we see an 
instantaneous measure of the number of stars in a particu- 
lar evolutionary phase (n), which should vary over time by 
~ \fn. 

The values in Table [4] are subject to biases in our de- 
tection statistics. These biases are dominated by three key 
factors: 

(i) The assumption of homogeneity across the core of the 
Sgr dSph, and that the observed GIRAFFE fields are rep- 
resentative of the core region. We see no obvious difference 
among the stars observed in the nine different VLT fields, 
either in the membership fraction of the individual regions 



in Table |4j nor the fraction of carbon stars in each fielcfl 
We note that M54, which comprises of a population of dif- 
ferent metallicity, lies within the core region, but this is rep- 
resents a small («4-5%) contribution to the total number 
of objects. We have not sampled M54 spectroscopically, so 
it should only have a secondary effect on the fractions we 
list in Table [4] but it does make a contribution to the total 
number of objects. 

(ii) The sampling of the GIRAFFE spectra within each 
field. While observations were limited by the placement of 
fibres, we have no reason to believe that the precise position 
within a field should bias our spectra. FLAMES targets were 
essentially chosen randomly from sources within regions (a)- 
(e) within the constraints of the GIRAFFE field. The UVES 
sample is slightly biased towards brighter targets, but these 
only represent a few percent of our targets and, because we 
can be confident of the spectral types and radial velocities of 
these stars, should not lend a significant bias to the results 
presented in Table [4] 

(iii) The difference in observation characteristics and 
spectral typing efficiency in the SAAO observations. As no 
radial velocities are measured for these stars, this can only 
bias the carbon star fraction. Within this, there should only 
be a significant impact on the fraction in region (b), as we 
assume all stars within region (c) are carbon stars and there 
are only eight stars successfullly observed by SAAO in re- 
gion (a). While our detection efficiency of carbon stars in 
region (b) may be biased due to the variable signal-to-noise 
of these observations, the small number of carbon stars this 
region can realistically host means any bias can have little 
impact on the total number of carbon stars in the core re- 
gion. Removing the SAAO spectra from the analysis leads 
to an increased error budget (primarily due to the absence 
of carbon stars in region (b)), but no substantial changes to 
the timescales involved, as the error budget is already quite 
large. 

We therefore do not expect biases in our data to affect the 
timescales we find by greater than a few percent and have 
included a conservative 5% error in our timescales to account 
for this. 

4-3.2 The stellar death rate 

The first step in our calculation is to determine the evolu- 
tionary rate near the top of th e RGB in kyr mag -1 . The 
Dartmouth evolutionary tracks (|Dotter et al.ll2008T ) contain 
a pre-computed model grid spaced at A[Fe/H] = 0.5 dex 
and AAf = 0.05 M Q . We assume solar [He/H]. Interpolat- 
ing over this grid, we identify the time taken for a star to 
traverse the upper-RGB region (f). Using bounds of -0.7 < 
[Fe/H] < -0.4 dex, 4 < age < 8 Gyr and assuming [a/Fe] = 
+0.2 dex, we find that stars should traverse region (f) over 
a period between 1020 and 1106 kyr. 

The next step is to determine the number of Sgr dSph 
stars that exist in this region. We focus on the galaxy's core 
(where our observations were taken) to determine numbers 
as our membership selection efficiency will vary across the 

5 We note that the overall fractions may vary from field to field 
due to the differing contamination by foreground stars, but this 
has little effect on our colour-selected samples. 
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Table 4. Star counts from the regions shown in Figure [T] 



Rej 


;ion Total count 


Confirmed 


Carbon star 




r.v. 


member? 


Implied number in 


core 




Body 1 


Core 1 


M-type Carbon 


fraction 2 


Yes 


No 


Fraction 2 


Carbon stars 3 


Members 3 


(a) 


7327 


2038 


182 2 


1.3 ± 0.7% 


234 


7 


97.0 ± 0.9% 


12-37 (12-38) 


1959-1995 


(1915-2039) 


(a) 4 


2385 


700 


119 2 


2.0 ± 1.0% 


114 


4 


96.4 ± 1.4% 


6-20 (5-20) 


665-685 


(640-710) 


(b) 


563 


162 


44 4 


9.1 ± 3.4% 


38 


3 


92.1 ± 3.5% 


4-14 (4-15) 


144-155 


(132-167) 


(c) 


161 


35 


14 


98.0 ± 2.6% 


4 





93.9 ± 8.2% 


35 5 (29-41) 


35 5 


(29-41) 


(e) 


17 


8 















0-8 (0-11) 


0-8 


(0-11) 


(f) 


993 


301 


57 1 


2.4 ± 1.6% 


55 


2 


96.0 ± 2.1% 


2-11 (2-12) 


283-295 


(267-312) 


(g) 


216 


56 


16 


2.5 ± 2.3% 


13 


1 


91.2 ± 6.1% 


3 (0 3) 


48-54 


(41-62) 



defined in Figure |3] 2 The error is defined using the probability and error approximations of Wilson ( 1927). 3 Thc left column of 
each pair denotes the implied number, the right column (bracketted) includes the Poisson noise necessary to accurately determine 
evolutionary rates. 4 Only stars above K s = 11.533 mag. 5 AU stars in this region are assumed to be carbon stars within the Sgr 
dSph. 



galaxy (see Figure [3] for the definition of areas). Table [4] 
shows we estimate that 267-312 Sgr dSph stars in region (f) 
of the galaxy's RGB. 

We must now remove the contaminating AGB from the 
RGB sample. We can estimate the number of these by count- 
ing the stars in the region above the RGB tip (region (g)), 
which should all be AGB stars. Table H lists 41-62 Sgr 
dSph AGB members in this region. Typically, models (in- 
cluding the aforementioned Dartmouth models) predict that 
the rate of evolution up the AGB increases with luminosity, 
however the observed 7\s-band luminosity function between 
9.58 < K s < 10.58 mag is flat (-11% ± 20% increase in num- 
ber per magnitude fainter in K B ). We can therefore expect 
a similar number of AGB stars to pollute the RGB sam- 
ple, albeit with an increased uncertainty due to the slope 
of the AGB luminosity function. We calculate that 32-65 
Sgr dSph AGB stars pollute the RGB sample, and there- 
fore expect that region (f) contains 202-280 Sgr dSph RGB 
stars. 

Based on the Dartmouth models' evolutionary rate 
through this period, we expect one star to evolve off the 
RGB every 3500-5700 years in the galaxy's core. We expect 
the majority of these stars to subsequently evolve through 
their horizontal branch and AGB evolution. Given it takes a 
zero-age horizontal branch star only ~ 100-200 Myr to com- 
plete its AGB evolution as far as the RGB tip (|Marigo et al.l 
2008), we can say that the rate of AGB evolution (in units 
of years per star) is similar to the RGB. Therefore, one star 
in the core 3° x 3° region also evolves off the AGB every 
3500-5700 years. Assuming evolutionary rates in the rest of 
the galaxy are similar (which is only true to first order due 
to the radially-changing age and metallicity of the popula- 
tion) then we can simply scale by the number of sources 
in regions (a-c) to find the evolutionary rate for the entire 
main body of the galaxy. Under these assumptions, a star 
would evolve off the AGB every 1000-1700 years somewhere 
in the main body of the galaxy (we remind the reader that 
the main body is defined in this context in Figure [3]). 



4-3.3 The thermal-pulsing and optically-obscured lifetimes 

We can further note that the upper AGB (above the RGB 
tip; regions (b) and (c)) contains «183 bona fide Sgr dSph 
AGB stars (core region only). Using the above timescale 
per star, we can estimate the time a star spends on the 
upper AGB to be between 630 and 1060 kyr, on average. 



Due to the spread in stars' initial properties, this is likely 
to vary quite significantly, probably by more than the error 
budget. The total extent of the upper AGB is AK S ~ 1.2 
mag. If we assume that AGB evolution is 3.26-6.83 x (viz. 
280-202 RGB stars / 41-62 AGB stars) faster than RGB 
evolution which, according to the Dartmouth models, takes 
1020-1106 kyr per 0.4 mag, this sets a duration to the upper- 
AGB lifetime of 430-1070 kyr. Similarly, if we assume that 
thermal pulses begin «0.3 7f a -band magnitudes below the 
RGB tip (in accordance with the above models) then the 
same evolutionary rate (in terms of kyr mag -1 ) gives a TP- 
AGB lifetime is 530-1330 kyr. 

Typical therm al pulse intervals for solar- mass AGB 
stars are ~100 kyr: iKarakas fc Lattanziol (|2007l ) determine 
143 kyr for a 1 Mq, [Z/H] = -0.68 model including mass 
loss, and 75 kyr for a 1 M , [Fe/H] = -0.38 which does 
not include mass loss. Taking a range of 75-150 kyr implies 
that these stars undergo 3-14 thermal pulses on the upper 
AGB, or 4- 17 thermal pulses on the TP- AGB in total. For 
comparison. iMarigo fc Girardil (|2007l ) predict a total of 11 
thermal pulses for a 1.0 M , [Z/H] = -0.38 model. 

Similarly, we can estimate the time a star spends in the 
optically-obscued carbon-rich phase by simply counting the 
stars in region (c). This region contains 35 stars, of which 
we believe all are carbon-rich members and all of which are 
sufficiently reddened in ( J — K s ) to be denoted as optically- 
obscured. Accounting for the Poisson error in this count, we 
find the average lifetime of this phase to be 100-235 kyr. 

These figures should be treated as approximations, how- 
ever, for two reasons: 

(i) We assume that net circumstellar absorption by dust is 
negligible in the K s band. The reddened, carbon-rich AGB 
extends slightly further above the AGB in its-band, and 
may extend somewhat further in bolometric luminosity, but 
contains relatively few stars, meaning stars do not typically 
survive long in this phase. This could potentially increase 
the maximum upper- AGB lifetime by another ~25%, based 
on relative star counts between regions (b) and (c). This 
would give a maximum of <17 thermal pulses on the upper 
AGB, or <21 thermal pulses on the total TP-AGB. 

(ii) We assume all stars within our colour cuts belong to 
the bulk population. Some will be part of the (very) metal- 
rich populations, though the number of such stars currently 
on the RGB and AGB is a factor of several smaller than 
the number in the bulk population. The metal-poor popu- 
lations are largely unaccounted for here: the giant branch is 
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Figure 9. Maximum C/O ratio attained by stars at a 
given metallicity and i nitial mass, based on isochrones from 
IMarigo fe Girardil l|2007l l. Light green areas show regions where 
C/O never exceeds unity; red areas show regions where stars be- 
come carbon rich, with darker colours representing increasing val- 
ues of the maximum C/O ratio. T he dashed lines show the ranges 
found in the Sgr dSph, based on ISiegel et al ] l|2007l l and using 
IMarigo fe Girardil l l2007ft to covert ages into masses. 



well-hidden within the foreground Galactic population. Only 
significantly-reddened stars (of which there will be very few) 
from this population will enter our colour cuts. This could 
change the ratio of stars in the different regions, hence our 
relative timescales, but only by a factor much less than the 
sub-populations contribution (by fraction) to the total pop- 
ulation. We therefore estimate this will have at most a few 
percent effect on our results. 

As our error budget is quite large, the net effects of these 
assumptions have little effect on our overall timescales. 

4.4 Which stars become carbon-rich? 

The origin of the Sgr dSph's carbon-star population has so 
far been unclear. Which of the populations do carbon stars 
come from? We have identified 23 carbon stars in our spec- 
tra (Tables [2] & [3]), all but four of which are new5 This 
more than doubles the number of confirmed carbon stars in 
the galaxy. The brightest identified carbon star , #2001, is 
the mass- losing star PQ Sgr (P=204 d, ASAS; IPoimanskil 
Il997h . As it was observed with SAAO, it has no radial ve- 
locity measurement. Its location on the colour-magnitude 
diagram, however, suggests that it is probably a Galactic 
star: it is the brightest carbon star in Figure [5] and sits three 
K a magnitudes above the carbon-rich AGB. 

Figure [9] shows the region in mass-metallicity space 
where stars are expected to bec ome carbon-rich, accord - 
ing to the theoretical isochrones of lMarigo fc Girardil (|2007l ) . 
This region is bounded at high mass by hot bottom burn- 
ing, which destroys ca rbon which would otherwise enhance 
the surface C/O ratio (IScalo et al.lll975l ; iDownes et al1ll996l: 
IVentura fc Marigol r2010). More importantly, it is bounded 
at low mass and high metallicity by the amount of dredge- 
up that occurs: low mass stars have little dredge-up, which 
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can prevent C/O reaching unity; high metallicity stars have 
high natal oxygen abundance, meaning more carbon needs 
dredged up to increase C/O to unity. 

All three populations of the Sgr dSph are near this low- 
mass-high-metallicity limit and models of this phase of evo- 
lution are not yet sufficiently accurate to say which side of 
this li mit each of the pop u lation s are o n. However, the mod- 
els of IMarigo fc Girardil (|2007f ) and iKarakas fc Lattanziol 
i|2007l ) both suggest the metal-intermediate population is 
least likely to produce carbon stars, while carbon stars may 
form in the metal-rich and (possibly) metal-poor popula- 
tions (note that Figure assum es a solar [q/Fe ] ratio) . This 
has already been suggested by lLagadec et ail ([2009), who 
claim the most-highly-reddened carbon stars are metal-rich 
on the basis of the C2H2 and SiC features in their infrared 
spectra. 

There are two important reasons why we now believe 
that the bulk population becomes carbon-rich, despite the 
above results. Firstly, we do not see any significant radial 
segregation of the carbon stars compared to the rest of the 
AGB population (Figure|3J, of which they represent a signif- 
icant fr action. Due to the strong metallicity gradient in the 
galaxy l|Siegel et al.ll2007h . such a radial segregation would 
occur if carbon stars only arise from a fraction of the popu- 
lation. 

Secondly, the stellar death rate correlates with the 
(carbon-rich) planetary nebula formation rate. Three asso- 
ciated planetary nebulae are known within the area we de- 
fine as the Sgr d Sph's main body: Hen 2-436, Wray 16-423 
and StWr 2-21 dZiilstra et all l2006al : iKniazev etafl 120081 ; 
lOtsuka et alj|201ll ). All have well-defined metallicities, and 
are known to belong to the metal-intermediate or metal-rich 
populations. Crucially, all are carbon-rich. The oldest, StWr 
2-21, has a dynamical age of 5200 years. This is consistent 
(within errors) with every AGB star in the bulk population 
becoming a carbon-rich planetary nebula, leaving a ~0.61 
M© white dwarf remnant. It is inconsistent with solely the 
much smaller, metal-rich populations producing planetary 
nebulae. 

On balance, it seems highly likely that at least the ma- 
jority of stars in the bulk population of the Sgr dSph become 
carbon rich and end their nuclear-burning lives by produc- 
ing a planetary nebula. We cannot currently comment on the 
eventual fate of the metal-rich and metal-poor populations. 



4.5 The S stars 

At low metallicity, less carbon needs dredged up to increase 
C/O to above unity. The S-star transition should therefore 
occur in a shorter period of time and thus few S stars are 
expected. However, the wavelengths of the dominant ZrO 
bands in these spectra are very close to those of the TiO 
bands (Figure [4}, making visual identification of S stars pos- 
sible only when the C/O ratio is very high (C/ O > 0.9, based 
on a c omparison to the models presented in Ivan Eck et al] 
(loll])). 

One of our S stars (#424; K s = 11.34 mag) is ^0.7 mag 
below the RGB tip (Figure [5}. This gives it a luminosity of 
~1000 L , placing it below the start of the TP- AGB, be- 
fore significant dredge-up of carbon-enhanced material has 
occurred. Two possibilities may explain its existence. It may 
be an extrinsic S star: a secondary star in a binary system 
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which has undergone mass transfer from a carbon-rich pri- 
mary, but not quite enough to make it carbon-rich. However, 
the location of a carbon star at a similar magnitude, and rel- 
ative lack of S or carbon stars at fainter magnitudes, suggests 
it has an intrinsic origin. If so, it is likely that both the faint 
S and faint C star are currently at a luminosity minimum, 
having recently gone through a thermal pulse. Carbon stars 
like these are seen in other environm ents, including the dis - 
rupted dwarf galaxy core, uj Centauri |van Loon et al.ll2007t ). 
While it is not clear whether the stars in to Centauri are in- 
trinsic or extrinsic either, they cluster near the RGB tip and 
are thus more likely to be intrinsic^. The remaining four S 
stars (#361, #909, #942 and #1054) are almost certainly 
intrinsic S stars. 

The brightest S star, (#361; K s = 9.28 mag) is above 
the obvious AGB tip (K s = 9.50 mag). Similarly to #424, 
there are three possibilities to explain this. Firstly, this may 
be a variable star caught in a luminous part of its pulsa- 
tion cycle (var i ations of AK S ~ 2 mag are not unusual; cf. 
iReikuba et ail (|2003r i). Secondly, the star may be currently 
undergoing a thermal pulse and thus be abnormally bright 
(though this phase is very short in comparison to the ther- 
mal pulse cycle). Thirdly, it may be a more-massive star 
in the metal-rich population, which survives to significantly 
higher luminosities on the AGB. 

If every AGB star in the bulk population does indeed 
become carbon rich, we can use the evolutionary rates in 
Section \4. 31 to determine the speed of carbon enrichment in 
that population. There are 700 objects brighter than the 
faintest S star (#424) and redward of our colour cut within 
the Sgr dSph core region. Of these, 114 are confirmed Sgr 
dSph members and 4 non-members. We can therefore ex- 
pect of order 5/114 x 114/118 x 700 « 30 (± 13) S stars in 
the Sgr dSph core. If one star evolves past a certain point 
every 3500-5700 years, the transition from C/O ~ 0.9 to 
C/O w 1.0 should take of order 60-250 kyr. This implies 
that a single dredge-up event (which happens every thermal 
pulse, or ~100 kyr) typically enriches the stellar surface with 
between ~4%-17% of the star's initial carbon abundance. 
This value may increase (slightly) further if star #424 is 
an extrinsic S star. Such values are not unexpected for l ow- 
mass, low-metallicity stars (e.g. Mari go fc Girardll2007l ). If 
the amount of carbon enrichment exceeds ~10% per thermal 
pulse, many stars may progress straight from the oxygen- 
rich M-star phase to the carbon-rich phase, without passing 
through the S star phase at all. However, if carbon enrich- 
ment per pulse is lower, stars might survive one or two ther- 
mal pulses in the S-star phase. 



4.6 The carbon-rich population 

4-6.1 Carbon richness and dust production 

With a single exception (#588, K B = 13.53 mag), all the 
carbon stars identified in the Sgr dSph are either within one 
magnitude of, or above, the RGB tip. This argues that the 
vast majority of carbon stars in the galaxy are intrinsic car- 
bon stars, excepting #588. The extent of the oxygen-rich 



AGB to ~1.2 magnitudes above the _K" s -band RGB tip fur- 
ther suggests that the transition to C/O > 1 typically (but 
not necessarily always) occurs above the RGB tip. Unsur- 
prisingly, this is also the region where thermal pulses are ac- 
tive, and the convective zone of the stellar atmosphere allows 
large amounts of carbon-enhanced material to be brought to 
the surface. 

The oxygen-rich AGB tip is defined by 2MASS 
18514859-3120157 at (J - K s ) = 1.30 ma g and K s = 9.50 
mag. Using the SED-fitting method of McDonal d et al 



(2009), incorporating the additions used in lMcDonald et al 



l|2011af ). we m odel the luminosity of this star using the BT- 
Settl models of lAllard et"aH (120031 ). assuming [Fe/H] = -0.55 
dex, M = 1.1 M , d = 26.3 kpc and E(B - V) = 0.14 mag. 
We arrive at a luminosity of 5200 L© for the AGB tip, with 
an associated error imparted by dust-enshrouding, stellar 
variability, uncertain distance and incomplete photometry 
estimated at ±500 Lq. With the exception of the aforemen- 
tioned S star (#361), all stars beyond this point appear to 
be highly-reddened carbon stars. 

Close to the RGB tip, stars also start to diverge from 
the main giant branch toward much redder (J~K S ) colours 
(region (c) of Figure QJ . This indicates the presence of dust 
production in these stars. We can note that all the observed 
stars which lie redward of the giant branch are carbon stars, 
which strongly suggests that all of the stars in region (c) are 
carbon starqj. This is perhaps unsurprising, as dust products 
of carbon stars (notably amorphous carbon) are much more 
opaque, especially in the near-infrared, than the dominant 
sources of opacity in oxygen-rich stars (silicates and possibly 
metallic iron). 

There also appears a population of dustless carbon stars 
lying within the main giant branch, but these stars are few in 
number compared to the dusty, reddened stars. It is possible 
that enhanced dust production begins to occur shortly after 
the star becomes carbon rich, though we cannot rule out 
that the carbon stars have dustless episodes (related, e.g., 
to their thermal pulse cycle). This pattern of a few dustless 
and many dusty carbon stars on the AGB is also seen in 
other (mostly d warf) galaxies with populations of similar age 
and metallicitv (IWhitelock et al.ll2009l ; iMenzies et"alll201ll ; 
iBover et al]|201ll ; Ijavadi et alj|201ll l 



4-6.2 Carbon star lifetimes 

We can compute the average time stars spend in the carbon- 
rich phase by performing a similar analysis to that in Section 
14.51 We note that individual carbon stars may have consid- 
erably different lifetimes than this average. 

We first require the number of intrinsic carbon stars in 
the Sgr dSph. As we discussed in the previous section, it is 
likely that all the reddened stars (defined by region (c) in 
Figure [TJ in the Sgr dSph are carbon stars. Whether region 
(e) also contains a carbon star population is uncertain, and 
we therefore have between 29-52 reddened carbon stars in 
the Sgr dSph core, as a temporal average. Again, we split 
the giant branch as before into regions (a) and (b). Above 



7 Note, however, that decreasing ZrO band strengt h with temper- 
ature makes hotter S stars more difficult to identify dZiilstra et al.l 
|2004) . 



8 Changes in molecular opacity in carbon stars also lead to red- 
dening compared to otherwise-identical oxygen-rich stars, but this 
accounts for, at most, a few tenths of magnitudes in (J — K s ). 
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the RGB tip, region (b) contains an estimated 4-15 carbon 
stars in this region (Table [3J. In region (a) (below the RGB 
tip), there are 700 2MASS objects brighter than our second- 
faintes10 carbon star (#1044, at K s = 11.533 mag). We 
estimate that there are 5-20 carbon stars in this region. 
Adding these contributions together makes between 38 and 
87 intrinsic carbon stars in the Sgr dSph core in total (note 
that, at the present epoch, we estimate there to be between 
45 and 77 such stars). 

If these are mostly intrinsic AGB carbon stars, and are 
representative of the bulk population, then our previous evo- 
lutionary rate (3500-5700 years) yields a typical carbon star 
lifetime of 130-500 kyr, or approximately two to five thermal 
pulses. 

We note that a carbon star lifetime of 130-500 kyr is 
much shorter than the typical TP-AGB lifetime (530-1330 
kyr; Section 14. 3|) . therefore its carbon stars should never 
reach C/O ratios much higher than unity (perhaps averaging 
C/O m ax ~ 1.2-1.5, on the basis of A C/O ~ 0.1 every ~100 
kyr (Section [43]). The Sgr dSph bulk population is evidently 
close to the limit for pr oducing carbon-rich stars (indeed, 
iMarigo fc GirardU |2007T ) predict the bulk population should 
not become carbon-rich; see Figure [9]). One may therefore 
find it unsurprising that carbon stars have a short lifetime, 
as they become carbon-rich shortly before the superwind 
would otherwise occur. 

However, it has been previously theorised (see Sec- 
tion [TJ that process of becoming carbon rich heralds the 
star's death. In this scenario, the increased opacity of car- 
bon molecules causes the star to expand, such that more of 
the extended atmosphere reaches the dust-producing zone. 
Meanwhile carbonaceous dust, being more opaque than 
silicaceous dust, provides increased dust-driving efficiency 
within the wind. This increases the mass-loss rate to such 
an extent that the star's atmosphere can be dissipated on 
a much shorter timescale. This scenario requires that the 
oxygen-rich progenitor star is already close to the threshold 
needed to generate a superwin d, which seems primar ily set 
by its pulsation characteristics (|van Loon et al.ll2008l ). This 
criterion may hold for the bulk population of the Sgr dSph, 
and it thus provides a useful laboratory in which to test such 
theories. 

4-6.3 Return of enriched material to the ISM 

We noted in Section 14.6.11 that all the stars with very red 
(J—K s ) colours appear to be carbon stars. This is unsurpris- 
ing when one considers that circumstellar carbonaceous dust 
provides a much higher (Aj-Aks) extinction than an equiv- 
alent amount of oxygen-rich dust, and that reprocessed emis- 
sion from this dust can already be significant at iG-band. 
Conversely, the increased transparency of oxygen-rich (sili- 
caceous) dust means that (Aj-Aks) will be low, even for sig- 
nificant amounts of dust. Oxygen-rich dust also tends not to 
re-emit effectively shortward of the Si-0 bending mode near 
10 /xm, unless it contains a significant fra c tion of iron (e.g. 
Kemper et all |2002| ; rVerhoelst et af]|2009l ; iMcDonald et al] 
20ld b 

9 Again, we assume that #588 is an extrinsic carbon star and 
thus ignore it. 
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Figure 10. 2MASS-WISE colour-colour diagram of the Sgr 
dSph. Black points show the entire 2MASS- WISE cross- 
correlated catalogue brighter than K B = 11.22 mag. Small, blue 
dots denote the subset of colour-selected AGB Sgr dSph stars; 
red dots denote observed carbon stars; green dots denote ob- 
served M-type stars. The lines mark the approximate locations 
of the carbon-dust-producing and silicate-dust-producing evolu- 
tionary tracks. The S stars, marked with magenta circles, all lie in 
the largely-dustless population near the confluence of these two 
tracks. 

Figure fTCTl shows a colour -colour diagram i ncorporating 
both the 2MASS and WISE (|Cutri et alj|2012f ) data. These 
colours have been chosen so that the heavily-reddened car- 
bon stars, which usually display very red colours at almost 
all near- and mid-infrared wavelengths, lie off to the right of 
the diagram, as indicated by the solid red lines. Meanwhile, 
the oxygen-rich giants have blue (J — K s ) colours due to the 
lack of dust opacity, but large ([4.6]-[11.6]) colours thanks 
to the 10-/xm silicate feature, as indicated by the dashed 
green lines. This diagram appears to separate the different 
stellar types reasonably well, but is not 100% effective. A 
full analysis of the relative mass loss from, and dust produc- 
tion by, oxygen- and carbon-rich stars is beyond the scope 
of this paper. However, Figure \W\ clearly shows that the car- 
bon stars are not the only dust producer in the Sgr dSph: 
the oxygen-rich, M-type stars return significant amounts of 
dust as well. 



5 CONCLUSIONS 

In this work, we have performed a spectroscopic study of the 
nearby dwarf galaxy, the Sgr dSph. We present 1142 spectra, 
estimating radial velocities for 1058 and finding with high 
(>99%) confidence that at least 592 are members. Highlights 
of these spectra include: 

(i) Spectra of 23 carbon stars, of which 19 are newly- 
determined; 

(ii) Spectra of 5 S stars, all of which are newly- 
determined; 

(iii) Spectra of ~2-6 stars which may belong to the metal- 
poor population. 

From these spectra and associated photometry we find: 

(i) The carbon- and oxygen-rich stars separate into two 
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distinct branches on colour-magnitude and colour-colour di- 
agrams, leading to the estimation that there are 31 ± 14 S 
stars and 63 ± 12 carbon stars in the central 3° x 3° of the 
Sgr dSph in total. 

(ii) The average star in this region is separated from its 
nearest evolutionary neighbour by between 3500 and 5700 
years of evolution. Extrapolating this to the portion of the 
galaxy not obscured by the Galactic Bulge gives one star 
per 1000-1700 years. 

(iii) We use the observed planetary nebula frequency 
within this region (3) to argue that the entire bulk pop- 
ulation appears to form carbon-rich stars, then planetary 
nebulae. 

(iv) On this basis, stars spend an average of 60-250 kyr 
as S stars (0.9 < C/O < 1) and 130-500 kyr as carbon stars 
(C/O > 1). 

(v) If every star becomes a carbon star, each thermal 
pulse should increase the C/O ratio by 0.04-0.17, and stars 
should typically reach a maximum C/O barely exceeding 
unity, though individual stars may reach higher C/O ratios. 

(vi) Carbon stars return a significant amount of carbon- 
rich dust to the ISM, though it is not clear whether or not 
this is enough to exceed the return from oxygen-rich stars. 

(vii) A small velocity gradient may exist, implying either 
tidal distortion or slow rotation around the galaxy's major 
axis. 

We suggest that further work is needed to determine el- 
emental abundances of these stars, with particular emphasis 
on identifying their metallicities. 
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